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Abstract 

Thermally activated escape of an over-damped particle from a metastable 
well under the action of a time-ramped force is studied. We express the mean 
first passage time (MFPT) as the solution to a partial differential equation, 
which we solve numerically for a model case. We discuss two approximations 
of the MFPT, one of which works remarkably well over a wide range of loading 
rates, while the second is easy to calculate and can provide a valuable first 
estimate. 

Introduction - Thermally activated escape of a particle from a metastable well, has 
found numerous applications in a variety of systems M. Escape under the additional action 
of a time-dependent force constitutes a non-trivial generalization. Several studies @-§J have 
been devoted to this problem for a sinusoidal force, which is particularly interesting since 
this system shows stochastic resonance 0. The topic of the present paper is the effect of a 
time-ramped force on the escape rate. Apart from its fundamental significance, motivation 
to study this problem arises from recent work on the dynamical strength of molecular bonds 
g-gg. The strength of a single bond can be measured experimentally using atomic force 
microscopy where one in practice applies a time-ramped force Pfl|TT|-|T3|j . Evans and co- 
workers HI! pointed out that the rupture strength of such a bond depends on the loading 
rate. This behavior has been seen not only in experiments but also in simulations of a 
Langevin equation ||,[10| and molecular dynamics simulations at very large loading rates 
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For the general problem of diffusive escape from a metastable well under a force that in- 
creases with time, we express the mean first passage time (MFPT) as the solution to a partial 
differential equation (pde). This exact approach differs from previous work that introduced 
an approximation based on instantaneous decay rates PJlQl. The numerical solution of the 
exact equation is then compared with both this and another simple approximation. The first 
one works very well over a large range of loading rates. The second one can be calculated 
easily and still gives a reasonable estimate which is never off more than 30% over the entire 
range of loading rates. 

MFPT as an exact solution to a pde - We consider the motion of an overdamped particle 
in a one-dimensional potential V(x) under the action of a time-dependent force f(t). The 
Langevin equation for this particle reads 
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d t x = -V'(x) + f(t) + C (1) 

where the stochastic noise obeys the usual correlations < ((t)((t') >= 25(t—t'). Throughout 
the paper we measure energy in units of ksT and time and length such that the diffusion 
coefficient becomes 1. The potential V(x) has a metastable well centered at 
saddle point at x s and an activation energy Q = V(x s ) — V(x m ). While our approach holds 
for any f(t), we will often specialize to a linear ramp 

fit) = fit (2) 

with loading rate fi. 

In the absence of the force (/=0), the MFPT T(x) that a particle originally at x needs 
to escape from this well obeys the equation |T5| 



&T(x) = {-V'd x + d 2 x )T(x) = -I. (3) 

Here, £t is the backward Fokker-Planck operator. The boundary conditions are T(xj) = 
T(x r ) = 0, where for a metastable well xi is far to the left of the minimum and x r far to the 
right of the saddle point x s . The explicit solution T(x) as obtained by simple quadratures 
of (^|) is well known |]T5] ? p~6|] . 

In the presence of a time-dependent force, relation @ cannot be used without modifica- 
tion. To obtain an exact equation for the MFPT, one can get rid of the time dependence of 
the right hand side in (JJ) by writing this Langevin equation as a system of two equations 

d t x = -V(x) + f(r) + C, 

d t r = 1. (4) 

In this form, the problem formally corresponds to a process homogeneous in time for the 
two variables x(t) and r(t). Therefore, one can apply the standard equation for the MFPT 
which becomes 

(d T + f(r)d x + ^)T(x,r) = -l. (5) 

The boundary conditions in x are as before T(xi, r) = T(x r , r) = 0. Since the deterministic 
variable r does not diffuse, the equation is first order in r. Therefore, one can specify 
only one boundary condition in r which we choose as T(x, r e ) = 0. For these boundary 
conditions, the solution T(x, r) gives the MFPT that the process (4) starting for t = at 
x and r needs to reach either one of the following boundaries: (i) the particle reaches the 
absorbing boundary at x\ or x r \ or (ii) the variable r becomes r e . If r e is large enough, the 
number of events in which the particle has not reached the absorbing boundaries at x\ or 
x r but rather that at r = r e are negligible compared to those in which it escapes over xi 
or x r . Therefore, in the limit r e — > oo, the solution T(x,t = 0) gives the MFPT for the 
ramped problem. The argument r = arises from the fact that the ramping starts with 
f(t = 0) = at t = 0. 

Escape rate for constant force - Before we present numerical solutions of (||), we discuss 
two approximations for the special case of a linearly ramped force f(t) = fit. To avoid a 
proliferation of symbols we use the same function symbol T for all MFPT quantities under a 
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time-varying force, it being unambiguous from the arguments which case is being discussed. 
Both approximations use the solution to the MFPT problem under constant force /, which 
we call T Q (x,f), to construct approximate solutions for the MFPT under the time-varying 
force fit) = [it. The MFPT T (x, f) obeys 

(fd x + tf)T (x, f)=-l, (6) 



which is easily solved by quadratures as [15 



Ux , /)=[/** r & f - r #r r r' r dy 



V(y) Jx ($/') J Xl ' Jx $>(y) J Xl <$>{y') . Ixi Jxi 

(7) 



with iff(x) = exp(— V(x) + fx). These integrals will be calculated numerically. 

A significant simplification occurs if one only requires the result ([^D for large force- 
dependent activation barriers 

Q(f) = Q - f(x s (f) - x m (f)) » 1. (8) 

A saddle point analysis of (|?D then yields the well-known Kramers rate 

^(/)=T (/)exp[Q(/)]. (9) 

The characteristic time 

T ° (/) = [V»(x m (f))V»(x a (f))]W (10) 

is the inverse attempt frequency. V" (x mj8 (f)) is the curvature of the potential at the mini- 
mum and the saddle point, respectively, whose locations depend on the force /. 

We use the numerical calculation of the MFPT under constant force (0), rather than 
the Kramers expression (|9|), to construct solutions for a time- varying force using two ap- 
proximations. Both approximations thus obtained are valid for all applied forces even when 
Q(f) » 1 no longer holds. 

Self- consistent constant force (SCCF) approximation - The idea of this crude approxima- 
tion is that the escape is dominated by a typical force which is determied self-consistently. In 
the solution (0) for constant force we write / = fiT (x, f) which leads to the self-consistent 
relation T = T (x, /J,T). Solving this equation yields the SCCF approximation T s (x, fi). The 
main virtue of this simple approximation is that it can be calculated easily. 

For large barriers, a further simplification can be achieved. Using the Kramers expression 
(|), the SCCF approximation T^(n) follows from solving the implicit equation 

Ts (M) = {V»(x m )\V»(x a W/* GXP(Q ° " (Xs ~ Xm)) (U) 

where we notationally suppress the / = [iTf dependence in x mjS = x m>s (f). 

Adiabatic approximation - A somewhat more refined approximation has been introduced 
previously within the context of stochastic resonance by Zhou et al || and for the ramped 
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problem in Refs. @,[T(J. While the SCCF approximation focusses on the typical relevant force 
for the escape, the adiabatic approximation incorporates history dependence. It assumes a 
time-dependent escape rate v(t). The probability that a given particle escapes at time t is 
then given by 

p(t) = v (t) expf- f dt'vit')]. (12) 
Jo 



Equation ( 12f) says that the probability of rupture is the product of two terms: the 
instantaneous rupture rate u(t) at a given time and the cumulative probability of survival 
up to that time. This approximation rests on the assumption that the escape process takes 
no time. In this limit it becomes exact. For the linear time-ramped process, one now replaces 
the time-dependence by a force-dependence via f = fit which leads to the relation 

p(f) = (1/aiM/) exp[-(l//i) f dfvtf')]. (13) 

Evans and Ritchie || discuss the conditions under which this probability has a peak at non- 
zero /*, which is defined as the rupture force. Since within our approach, only the MFPT 
is available, we will discuss this quantity 

roc rf 

T a {x,n) = {l/n) #exp-[(l//i) / dfvtf')]. (14) 
Jo Jo 

Since u(f) = 1/T (x, /) with Tq(x, f) from (0), the MFPT T a (x, //) in the adiabatic approx- 
imation can be obtained by performing the two integrals numerically. 

Numerical solution - In this section, we present numerical data which compare the MFPT 
obtained from the two approximations with the numerical solution of the exact equation (^j). 

Specifically, we use the potential 

V(x) = (27/A)Q x 2 (l-x). (15) 

This potential exhibits a metastable minimum at x m = with a saddle point at x s = 2/3 
and an activation energy Q . 

Fig. 1 shows the MFPT for two different values of the activation energy as a function of 
loading rate fi as calculated from (|5|). For fi = 0, the MFPT approaches the zero-force value 
as obtained from (|3|). The inset shows that over a large range T ~ u~ a where a ~ 0.8 is an 
effective exponent masking an extended crossover from a = 1 (with logarithmic corrections) 
to a = 1/2 as discussed below. 

The quality of the two approximation is shown in Fig. 2, where we present the ratio of 
the MFPT 

R a , s (x, (jl) = T ajS (x, (jl)/T(x, fi) (16) 

as calculated from either approximation to the exact value. 

The gross features of these data can be understood by identifying three regimes P,|i~0 . 
For small /i, one expects naively that the MFPT deviates from the value at zero force if 
H > Ji where 
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Jjl = Q(0) exp[-Q(0)]/[r(0)(x a - x m )}. (17) 

While this (exponentially small!) loading rate also signifies the range where the SCCF 
starts to deviate from the exact data, the adiabatic approximation holds well even for larger 
loading rates. This is remarkable since as shown in the Appendix systematic perturbation 
theory reveals that both approximations fail to reproduce the 0(/x) correction to the zero 
force result for a general potential. It seems that only in the additional limit of a large 
activation energy, Q(0) >> 1, is good agreement between the exact result and the adiabatic 
approximation restored. A mathematical proof of this feature is beyond the scope of this 
paper. 

For large fi, both the potential and the diffusion can be ignored and then the ratio 
displayed in Fig. 2 can be calculated analytically. Even though the process in this regime is 
then no longer thermally-activated, we include it for two reasons. First, the MFPT is still 
mathematically well-defined and can even be calculated analytically in the limit \i — > oo. 
Second, for the bond rupture problem described in Refs. |J and JTIJ, this regime has some 
interest, where it has been called "ultrafast". The time-ramped equation of motion becomes 
d t x = fit with the solution x(t) = x + fit 2 /2. This leads to the exact asymptotic result 



T(x,/i) « y/2( Xr -x)/n (18) 



and hence to a = 1/2 as mentioned above [ 17fl . 



For the SCCF approximation, the equation of motion becomes dtx = /, which leads to 
T(x) = (x r —x)/f. Replacing / by fiT(x) yields the relation 



T 8 (x,fi) ~ \J [x r x)/fu (19) 

and hence R s (x, //) ~ v^2/2 for large \x. In the same limit, the adiabatic approximation can 
be calculated by performing the integrals in fl"H]). One obtains R a (x,fj,) ~ \^n/2 which is 
closer to the exact result than the SCCF approximation. 

In an intermediate fi— regime, the behavior can be extracted most easily from the SCCF 
approximation. From the relations lnT,f = lnT x (0) — fiT^(x s — x m ) one obtains ~ 
ln/^//x and thus a = 1 with logarithmic corrections. Expressed in terms of a mean rupture 
force /* = [i, this behavior corresponds to /* ~ In // |J . As Fig. 1 shows, the data show 
a large crossover with an effective exponent a ~ 0.8. 

Conclusions - In this paper, we have analysed diffusive escape over a barrier in the pres- 
ence of a time-dependent force. We have derived an exact pde for the mean first passage 
time and solved it both by numerical integration and by use of two approximations. Com- 
parison of the (numerically) exact solution with the adiabatic approximation introduced 
previously shows that this approximation holds remarkably well for a large range of loading 
rates. Formally, our exact equation for the MFPT can easily be generalized to include space- 
dependent diffusion coefficients. Likewise, both motion in higher dimensions and inclusion 
of inertia terms (for smaller friction) are amenable to the same treatment. The numerical 
solution of the corresponding equivalent of (g), however, will become quite time-consuming. 
It will then be helpful to explore the two approximations whose virtues we discussed for our 
model case. 
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APPENDIX A: SMALL LOADING RATES 

For small //, perturbation theory of the SCCF approximation works out as follows. We 
first expand the solution of @ as 

T (x, f) = T®(x) + /T«(x) + 0(f). (Al) 

The first order term obeys &T^'(x) = —d x T^°'(x). We write its solution with the absorbing 
boundary conditions as T^(x) = A[—d x T^}. The linear operator A can be expressed easily 
by quadratures and is formally the inverse to . The final step consists in replacing / by 
uT^(x) = fj,A[—l). The first order result of the SCCF approximation thus reads 

T s (x,fi) = T(°>(x) +fiA[-l] A[-d x T^]+0(u. 2 ). (A2) 

A similar perturbative expansion for the adiabatic approximation shows that the term lin- 
ear in fi coincides with the SCCF approximation whereas the terms 0(fi 2 ) differ in both 
approximations from each other. 

Time (or rather r)-dependent perturbation theory in fi for the solution of the exact 
equation fl5|) can be set up as follows. Inserting the ansatz 

T(x,t;») = T®(x) +/i[T (1) (a;) + rT 1 (1) (x)] + 0(/i 2 ) (A3) 

in (|), yields in 0(/i, r 1 ) the equation C^T^ = —d x T^(x) with the solution T± = 
A[—d x T^]. Inserting this solution into the equation 0(fi, r°) yields C^T^ = — (x) with 
the solution Tq(x) = A[—A[—d x T^]]. Thus, the full equation leads to the perturbative 
result 

T(x,fi) = T®(x) - iiA[A[-d x T®]]. (A4) 

Since ( |A4| ) differs from (|A^), neither approximation, perhaps somewhat surprisingly, repro- 
duces the correct amplitude for the term linear in fi. Note, however, that the u. = value 
is reproduced by both approximations. We suspect that by considering the limit Q — > oo 
even the linear term shown in (|A^) approaches the exact one shown in (|A4]), but we have 
not yet been able to show this mathematically. 
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FIGURES 




-4 




ln( n) 



FIG. 1. Mean first passage time T(0, fx) as a function of loading rate [i for Q = 7 and Q = 5 
in the potential (|l5|). Inset shows the same data logarithmically. For all [i, the starting point has 
been chosen as x = 0. The absorbing boundary is x r = 1.6 and the (quite irrelevant) left boundary 
is at x\ = —1.6. 
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FIG. 2. Mean first passage time T(0, /x) as a function of loading rate /i for Q = 7 and Q = 5 
in the potential (|l5[). Inset shows the same data logarithmically. For all /i, the starting point has 
been chosen as x = 0. The absorbing boundary is x r = 1.6 and the (quite irrelevant) left boundary 
is at = —1.6. 
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